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Abstract: Magnetic field generated by the pulse-current flowing into air- 
core reactor of the experimental advanced superconducting Tokamak (EAST) fast 
control power supply may cause a thermal problem in output cabinet. To solve the 
thermal problem, a layer or two layers of shielding are arranged around the air core 
reactor. When the maximum amplitude of the pulse current is 1 500A, the calculation 
method of magnetic field distribution around the air-core reactor is studied, as well 
as the output cabinet of eddy current loss. Taking eddy current loss as a main thermal 
source, three dimensional (3D) magnetic-thermal coupled finite element models of the 
reactor are established according to electromagnetism and thermodynamics theories. 
Temperature field of air-core reactor is analyzed and discussed in this paper. The 
results show that the magnetic field around air-core reactor and the temperature of 
output cabinet are decreased under shielded condition. 
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1 Introduction 


Experimental advanced superconducting 
Tokamak (EAST) fast control power supply is a 
large capacity single phase inverter power supply". 
The system consists of six H-bridges branches, and 
each branch has an air-core reactor. In this system, 
the reactor is employed to limit current, steady 
voltage and compensate reactive power. Recently, 
the analytical method has been used to calculate and 
analyze the magnetic field distribution in air-core 
reactor without shield’*!, Electromagnetic shielding 
is a technique that reduces or prevents coupling 
of undesired radiated electromagnetic energy into 
equipment, so as to keep its operating compatibly 
with its electromagnetic environment. In order to 
shield the magnetic field, many scholars study a wide 
variety of shields, shields of different shapes and the 
number of layers'**!, In practical applications, air-core 
reactor during operation is prone to faults, including 
partial discharge, overheating and burnout'’*!, The 
temperature rise is one of the main aging mechanisms 
that have to be considered, especially for the overall 
system; therefore, installing a shield for reactor is 
attracting more and more attention", 

With regard to shield design problems, optimization 
is usually minimization of the field in a predefined 
area — the target region — and minimization of 
electromagnetic losses in a shield. However, it is 
also useful to include thermal aspects: a shield may 
obtain a high temperature because of electromagnetic 
losses and heat radiation from the heat source!'”"*), 
Even if the average temperature rise satisfies design 
requirement, the hottest-spot temperature in reactor 
might exceed the maximum temperature limits of 
insulation materials. Moreover, thermal field was 
simulated by finite-element method (FEM) after 
calculating the heat transfer coefficient by using 
nusselt number'*'*!, In calculation, the fluid dynamic 
behavior cannot be accurately described without using 
simulating fluid field. The calculation results were 
very different from the measurement results, so the 
heat transfer coefficient must be repeatedly modified. 


Therefore, a coupled analysis of magnetic field and 
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thermal field is mandatory to compute the temperature 
rise in the design stage. 

The purpose of this paper is to present a 
magnetic-thermal coupled model to accurately and 
efficiently predict the temperature in output cabinet, 
and use special shield made of composite materials 
to absorb strong magnetic field generated by the 
reactor. Different shielding materials will bring 
different results. When single-shielded, the effect of 
DT4A is the best and the effect of DT4A/GO is the 
best in two layers of shield. Temperature in single 
and double shielded condition reduces 38.1% and 
64.2%, respectively. Comparing the simulation and 
measurement temperature, it can be seen that the error 
is within the controllable range. Therefore, thermal 


problem of the cabinet can be solved in this way. 
2 Formulation 


2.1 Structure of air-core reaction 


The air-core reactor is mainly used for power 
supply rectifier which is intended to resolve the 
problem of inconsistencies in parallel rectifying 
element pressure drop. The reactor is made of 32 turns 
of winding copper coils and each coil is encapsulated 
glass fibers impregnated with epoxy resin. 

This structure enhances electric strength, 
insulation strength and mechanical strength of the 
air-core reactor. Simulation and actual models of 
the reactor are shown in Fig.1. Three reactors in the 
middle of the cabinet, the distance between the cabinet 
and the reactor is 240mm and the distance between the 


adjacent reactor distances is 330mm. 


Air-core * 


Shield 


Fig.l Models of the air-core reactor 


2.2 Analytic calculation of magnetic field around 


air-core reactor 


Fig.2 shows the irregular shape of each copper 
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coil, it is approximated as rectangular and composed 
of four elliptical arc edges. So the magnetic field of 
air-core reactor can be obtained by analyzing the 


magnetic field of each coil. 


Idl=Idxe, T> (a, -b) 


Fig.2 The model of single copper coil 

For calculating magnetic field around air-core 
reactor, its motion equation is given by Biot-savart 
law. All parameters are shown in Fig.2. 

A pe (1) 
An? iy? 
Where uo denotes the permeability of free space, /d/ is 
current element. 

(1) Magnetic field of rectangle line current. In 
this paper, the axis is the longitudinal direction of the 
magnet, and the current is defined as positive when it 
flows towards the positive x axis. It should be stated 
that for the rectangular winding the magnetic flux 
density will be composed from the shares of each of 
the rectangle side. It is assumed that the length d/ of 
the current element is infinite small. The unit vector 
between point S(x, -b, z) on the wire 1 and arbitrary 
point P(x, Y» z;) can be expressed as 


adla (x, -x)e, +(y; +b)e, +(z;-z)e, 2) 


Sr J@n +O; +b) +l -2)° 


Where e,, e, and e, denote the unit vector along the 
x, y and z axis. 
According to Eq.(1) and Eq.(2), the induced 


magnetic field at arbitrary field point P is expressed as 


_ My Idxe, xe, 
k 4n ph 
E Uoldxe, x[(x; -x)e, +(y; +b)e, +(z;-z)e,] 
O Anil -x Ob l-z 7? 6) 


ChineXiva ERAF 


So the x, y,» z; component of the magnetic flux 


density is 
dB,,, =0 
Ul (Zz; -Z) 
ly; = ° j? dx (4) 
~ Ani(x, =x) +(y, +b) +(z,- z)’] 
Mol (y: +b) 


dB, = pe dx 
 4nilx -x) +(y, +bY +z- ] 
The total magnetic field at point P generated by 


wire | is given as 


2 af Mol (z; = 2) dx 
w “nila -x +y +b +(2,-2) T 


lI (z, -z) — — 


an[(y, +b) +(2,-2) IL l 6) 
B, =f" Hol (yi +b) — dx 
í “Anl(x, —x) +y +b) +z -zT 
o mob aa, era 
An[(y; +by +(Z, -z)'] n i 
With 
= (x, -a} +; +b) + -2 
(6) 


n = (x, tay +y, +b) + (z-z) 
Similarly, magnetic field at point P generated by 
wire 2, 3 and 4 can also be calculated. So the magnetic 


field of rectangle line current is 


mia-  [b-») +y], 


j An (x; -ay +(z; -z)’] B n 
ma-a foD erv 
An[(x, +a) +(z,-zyY]| 4% 
Mol (z, - 2) (x-a) (x; +h, 
A 4ni (y, +b)’ +(z; -z)’] n 


mI(z; -z) 2. x), = 
An[(y,-b)? +(z,-2)"] 


uyl(y, +b) 2. x) m 


~ Anf(y, -b} +z, -2)] n =), 


B 


(7) 
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With 


r= f(x, =a)? +(y, =b)? + (Z,-2)° 
r= (x, +a)" +9, -bY +z, -2)" 


(2) Magnetic field of oval line current. Fig.2 


(8) 


shows an oval current-carrying coil, and the oval 


equation is 
x? y 
+-—=] 9 
(a+cy b? ( ) 


Through the transformation, solving the magnetic 


field distribution of elliptical coils can be converted 
into solving circular coil magnetic field distribution. 
Assuming elliptical shape factor k= b/a and X=k'"x, Y= 
k?y, Eq.(9) becomes 


X’ +Y’ =(a+c)b (10) 
With 
R=(a+c)b (11) 
Some equations can be obtained by Fig.2. 
r, = 1, (sin Oe, + cos Oe, ) 
R= R(cos Ge, +sin Ge, ) 
(12) 
r° =r} +R -27, sin Ocos o 
di = R(-sin ge, +cos ge, )dp 
In actual space, the magnetic flux density of any 
point can be obtained by the series of limited items. 
According to Eq.(1), Eq.(12) and complete elliptic 


integrals, the magnetic field of oval line current is 


given as 
_ MolRn Li 2Rr, sin 0 
On; (r +R?” a A +R? 
mIR MIR Sp, 2Rr,sin 0 \” 
A +R)? URRY A r +R 
u IRT, cos 0.) 2R sino y 
In +R’ A i +R? 
(13) 
With 
Y =/x; +(z,-z) 
© 1  4n=1(4n=3)+-3 (2n-1)(2n-3)-1 a 
" (2n-1)! gee) 2n(2n-2)--2 
4n-1)--3 (2n-1)(2n-3)---1 
D, = 1 (4n+1)(4n-1)---3 (2n-1)(2n-3) EN 
2n g” 2n(2n-2)--2 
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With the analysis above, magnetic field around 


air-core reactor is 


2 
ares, +B. +B? 4 
32 T ar : a 
Bey . (15) 
i=l = Toly +B. 
T i 
With 
a, = arccos ga ae 
' Va? +b? fa? + (bc) ie 
ee a +(a+ce) +2(b-c)y -2 
2a’? +(b-c)y (a+c +(b-cy 
Where B, B,., is magnetic flux density in the 


B, and 


ry; 


direction of each axis of elliptical coil, B, 
B 


12; 


is magnetic flux density representing the direction 
of the rectangular coil. 

For Eq.(15), the flux density at any point in 
actual space can be approximated by the series 
of limited entry and the magnetic flux density 
distribution of output cabinet can be calculated. So the 
eddy current loss in output cabinet can be analyzed 
by electromagnetic induction law and Ohms law. The 


equation is expressed as 
T? 
P, = —(JBY Arelh (17) 
6p 


Where Ape, h 


length of the output cabinet respectively, p is iron 


, and / are thickness, height and 
conductivity. 

(3) Strong coupling of magnetic-thermal field 
calculation. Fig.3 illustrates the Joule heat taken away 


from output cabinet by the cooling modes including 


Output cabinet 


Natural convection 


Thermal radiation 


Fig.3 Axisymmetric model of air-core reactor 
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heat conduction, natural convection and thermal 
radiation. The air-core reactor consists of 32-layers 
copper coils, some of which are shown in Fig.3. From 
the quantitative analysis point of view, it is necessary 
to establish a mathematical model in which the 
magnetic field equations, the heat conduction equation 
and the Navier-Stokes equations are linked through 
coupled variables. 

For the system of the reactor shown in Fig.3, 
there are three heat dissipation modes, involving heat 
conduction, natural convection and thermal radiation. 
For the output cabinet, the main heat dissipation mode 
is heat conduction. For the surface between output 
cabinet and surrounding air, the main heat dissipation 
mode is natural convection and thermal radiation. 
Finally, the temperature distributions of the output 
cabinet should be concerned. 

1) Heat conduction. For the output cabinet, the 


steady state heat conduction equation for solid is 


ð 422), 2 aTh wot +0=0 (18) 
ox\ ox} dy\ ody) dz\ dz 


Where A is the coefficient of heat conductivity, Q is 


given as 


the heat generation of unit volume in iron conductors 

2) Natural convection. The natural convection 
of air satisfies the Navier-Stokes equations, which 
consist of three groups of equations. For three 
dimensional incompressible steady fluids, the Navier- 
Stokes equations in Cartesian coordinate system can 
be simplified as following 


Continuity equation 


A(pu) | (pv) pw) _ 4 (19) 
ax ay oz 


Momentum conservation equations 


Ox oy 

ðw w > OP 

u— + pv— + pw— = UV w-—+ 
i ðx P ä 0 Oz f 


Energy equation is given as 


=AVT +O (21) 


u— +v—+w— 
Ox oy Oz 


e2 ƏT ƏT 
pe 
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Where p is fluid density, u is viscosity coefficient, c is 
the specific heat, u, v and w are the respective x, y 
and z direction potentials of the fluid speed, P is the 
fluid pressure, T is the fluid temperature, f,, f, and f, 
are the sources in the respective x, y and z directions. 

3) Thermal radiation. For a system of two surfaces 
(surface i and j) radiating to each other, the heat transfer 


rate between surfaces i and j is expressed as 
Q; = 0€,B,,A, (T$ a T) (22) 


Where ø is Stefan-Boltzmann constant, g, is the 
effective emissivity of surface i, p, represents the 
radiation view factor between surface i and j, A, is 
the area of surface i, and 7,, T, represent the absolute 


temperature of surface i and j. 


3 Calculations of a air-core reactor 


(1) Analysis of shielding thickness. To limit 
the leakage magnetic field inside the shield and use 
materials efficiently, the thickness of the shield should 
be considered. Penetration depth of the shield is 


[ 1 
ô= we (23) 


Where u is permeability of conductor, ø is 


expressed as 


conductivity of conductor. 

Take a reference line through the shield, it can 
be seen that the depth of the shield is based on the 
magnetic density of the reference line. Where the 
magnetic flux is high, the magnetic field passes. This 
distance is the penetration depth of high magnetic 
flux. Fig.4 shows penetration depth and magnetic 
flux density of the shield when AC amplitude and 


Magnetic flux density/mT 


A range of flux density 


0 10 20 30 40 50 60 70 
Distance/mm 


Fig.4 Penetration depth of the shield 
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frequency are 1 500A and 100Hz, respectively. 
Obviously, the penetration depth of the magnetic 
field shield does not exceed 20mm. Eq.(23) shows 
that penetration depth increases as the frequency 
decreases. To ensure that the magnetic field does not 
penetrate the shield, thus thickness of the shield is 
selected as 20mm. 

A shield is a metallic partition placed between 
two regions of space. It is used to control the 
propagation of electromagnetic fields from one region 
to the other. Shield changes not only the size of the 
magnetic flux density but also the direction of the 
magnetic field. Fig.5a shows the magnetic flux density 
distribution in no-shielded condition. It can be seen 
that the magnetic field around the four corners of the 
reactor is the strongest, and the center is the weakest. 
Because magnetic line is a loop formed from the 
inside to the outside. The magnetic flux densities 
around output cabinet in single shielded and double 
shielded condition are shown in Fig.5b and Fig.5c, 
respectively. Obviously, magnetic flux density is 
reduced in shielded condition and this phenomenon 
is more obvious when the double layer metal plate 


shields. When using the double shielding metal plate 


Flux density/mT 


l >~ 600 
400 
200 pgo 


(a) Distribution of magnetic flux density in no-shielded condition 


Flux density/uT 


600 


(b) Distribution of magnetic flux density in single shielded (DT4A) condition 
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(c) Distribution of magnetic flux density in double shielded 


(DT4A/GO) condition 


Fig.5 The magnetic flux density distribution 


(DT4A/GO, DT4A is closer to the field source), the 
shielding effect is better. 

(2) The influence of shield. Fig.6 represents the 
schematic drawing of simulation setup for simulating 
an air-core reactor shielding system. The reactor is 
surrounded by a shield, and the shield is between the 
reactor and output cabinet. The reference point is 
taken from the coordinate origin and multiple sets of 
data can be obtained from the reference point to judge 


the performance of shield. 


Air-core 
reactor 


Output cabinet 


Fig.6 Schematic drawing of the simulation setup 


The shield can be in any position between the 
reactor and output cabinet. So, the distance between 
shield and reference point can be changed. 

Magnetic field was measured from the reference 
point on the output cabinet. Fig.1 shows the linear 
relation between cabinet and the resultant magnetic 
field with shields. The shielding materials were grain- 
oriented electrical steel (GO), electrical pure iron 
(DT4A) and 1008 steel. All the materials are 20mm 
thick, and take amplitude (1 500A) of current which 
flow through air-core reactor. As a surrogate of 


shielding performance, shielding factor (SF) is used. 
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It is defined as the ratio of shielded magnetic field to 


unshielded magnetic field as follows 


B unshielded 


SF = 201g (24) 


shielded 


The greater the value of SF, the better the 
shielding effect. 

For magnetic fields generated by the reactor, 
the variation of the shielding factor as a function of 
distance is shown in Fig.7a. The reference point is 
fixed, and the metal plate moves toward the z axis 
at every 20mm distance. A noticeable reduction of 
magnetic fields by the shields is observed when the 
current (1 500A, 30kHz) flows through the reactor. 
When installing GO, matching the direction of 
grain orientation with field direction is important. 
Longitudinal direction of the GO sheets is an easy 
magnetization direction, whereas DT4A and 1008 
steel are magnetically isotropic. In the case of single 
layer shielding, single layer shielding factor (SF,) of 


the three magnetic materials decreases with increasing 


90 T T i T T T T T T T T T T T T T 
gol ^ —=—DT4A 

t —e—GO 
70 —A~ 1008 steel 


0 20 40 60 80 100 120 140 160 180 200 
Distance from the reference point/mm 


(a) Relationship between SF, and distance 


—sa— DT4A/GO —e—GO/DT4A 
—4&— DT4A/1008 steel—w—GO/1008 steel 
—- 1008 steel/DT4A—< — 1008 steel/GO 


SF,/dB 


L 1 1 L 4 L L 2 L 4 L 4 L 4 L 4 
0 20 40 60 80 100 120 140 160 180 200 
Distance from the reference point/mm 


(b) Relationship between SF, and distance 


Fig.7 Relationship between SF,, SF, and distance 
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distance. In the three materials, shielding effect of 
DT4A is the best and the maximum SF, is 105dB. 
Whereas, the shielding effect of 1008 steel is the best 
and the maximum SF, is 54.3dB. Shielding factor of 
DT4A compared with 1008 steel is higher than 48.2%. 
Taking into account all aspects, DT4A shows the most 
effective shielding performance at the electric current. 
Fig.7b shows the relationship between SF, (double 
layer shielding factor) and distance. When distance 
between the double layer shield and the reference 
point is not the same, the size of SF, will be changed. 
The distance from the reference point of shield can 
be moved as needed, as shown in Fig.7b. SF, is 
inversely proportional to the distance, the smaller the 
distance, the greater the SF,. Although the six kinds of 
results are similar, but the shielding effect of DT4A/ 
GO (DT4A is closer to the field source) is better. 
Therefore DT4A/GO as a shielding material is more 
appropriate. 

EAST fast control power provides a current and 
the range of current frequency is 100Hz to 30 kHz. 
Fig.8a and Fig.8b show single and double layer 
shielding factor of different frequency respectively. 
In this case, it can be seen that no matter what kind 
of combination, with the increase in frequency, 
SF, and SF, are gradually increased. Among these 
groups, DT4A/GO has the best shielding effect. 
Maximum shielding factor can be reached 121.45dB. 
1008 steel of the performance is the worst and the 
maximum shielding factor is 72dB. With the change of 
frequency, the double layer shielding method is more 
effective for shielding electromagnetic field. 

This result demonstrates that closely adjacent 
layers act as one body under the regime of shunt 
mechanism of magnetic shielding. For example, SF 
of closely adjacent DT4A/GO is simply dependent on 
an averaged permeability of DT4A and GO. On the 
other hand, the DT4A effectively reduces the strong 
magnetic field first, and then the GO sheet effectively 
shields the weakened field. In the case of 1008 steel/ 
GO pair, shielding performance was not improved 
since a strong magnetic field near the source was too 
high to be effectively shielded by 1008 steel. From 
the above analysis, it can be concluded that DT4A is 
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(a) Single layer shielding factor of different frequency 
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(b) Double layer shielding factor of different frequency 
Fig.8 Single and double layer shielding factor of 


different frequency 


the best one in terms of distance factor and frequency 
factor. 

(3) The temperature and eddy current loss of 
output cabinet. The temperature and eddy current 
loss of output cabinet can be simulated by ANSOFT 
and ANSYS Workbench. First, the eddy current 
density of the cabinet can be calculated. Then, the 
results will be imported into ANSYS Workbench. 
Finally, temperature distribution can be obtained by 
the two software and adding material. In no-shielded 
condition, the eddy current loss and temperature are 1 
225.2W and 68.5'C. Tab. shows the eddy current loss 
and temperature of the output cabinet under different 
shielding materials when pulse current flowing into 
air-core reactor. It can be seen that temperature in 
shielded condition reduces by 64.2% of no-shielded 
condition. 

The temperatures of the cabinet in shielded 
(DT4A/GO) and no-shielded condition are compared 
in Fig.9a and Fig.9b. As can be seen from the 
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Tab. The eddy current loss and temperature of the 


output cabinet under different shielding materials 


Material Eddy current loss/W Temperature/‘C 
GO 922.9 51.6 
DT4A 758.4 42.4 
1008 steel 1 087.5 60.8 
DT4A/GO 438.2 24.5 
GO/DT4A 499.0 27.9 
DT4A/1008 steel 550.9 30.8 
1008 steel/ DT4A 613.5 34.3 
GO/1008 steel 581.3 325 
1008 steel/GO 620.6 37.6 


temperature distribution diagram, the four sides 
of the temperature are relatively high, and this is 
exactly the case with the distribution of flux density. 
Simulation of the temperature in no-shielded condition 
and shielded condition error are 4.5% and 2.2%. 
Comparing the measured and simulated data, it can 
be seen that the thermal problem can be effectively 


solved by using shield to absorb strong magnetic field. 


Temperature/°C 
65.418 Max 
97 


Min 


(a) Distribution of temperature in no-shielded condition 


Temperature/‘C 
23.961 Max 
23.923 


23.661 
23.624 Min 


(b) Distribution of temperature in shielded (DT4A/GO) condition 
Fig.9 Distribution of temperature in no-shielded and 


shielded condition 


5 Conclusion 


In this paper the magnetic field and thermal 
field by experiment and FEM are analyzed. Magnetic 
shielding performances of three shielding materials 
in the vicinity of a center line on the output cabinet 


have been evaluated. In shielding with a single layer, 
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DT4A shows the best shielding performance. DT4A/ 
GO is the most effective shield in shielding with 


two layers. In sum, the minimum temperature of the 


cabinet dropped to 24.5'C by this measure. Finally, the 


thermal problem is effectively solved. 
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